Three-body forces and proton-rich nuclei 
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We present the first study of chiral three-nucleon (3N) forces for proton-rich nuclei along the N = 8 
and N = 20 isotones. Our results for the ground-state energies and proton separation energies are 
in very good agreement with experiment where available, and with the empirical isobaric multiplet 
mass equation. A future measurement of the ground-state energy of 22 Si will provide an especially 
sensitive test of chiral 3N forces. We predict the spectra for all TV = 8 and TV = 20 isotones to 
the proton dripline, which agree well with experiment for 18 Ne, 19 Na, 20 Mg and 42 Ti. In all other 
cases, we provide first predictions based on nuclear forces. Our results are also very promising for 
studying isospin symmetry breaking in medium-mass nuclei based on chiral effective field theory. 
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Exotic nuclei with extreme ratios of neutrons to pro- 
tons can become increasingly sensitive to new aspects 
of nuclear forces. This has been shown in shell model 
studies with three-body forces for the neutron-rich oxy- 
gen [l|, [3| and calcium [§] isotopes, which present key re- 
gions for exploring the evolution to the neutron dripline 
and for understanding the formation of shell structure. 
Three-nucleon forces predicted an increase in binding of 
the neutron-rich 51;52 Ca isotopes compared to existing 
experimental values, which was recently confirmed by 
high-precision Penning-trap mass measurements [J]. In 
addition, in this region, there are very promising devel- 
opments in coupled-cluster theory Q . 

Proton-rich nuclei provide complementary insights to 
strong interactions, exhibit new forms of radioactivity, 
and are key for nucleosynthesis processes in astrophysics, 
such as the rapid-proton-capture process that powers 
X-ray bursts 0, Q. Although the proton dripline is 
significantly better constrained experimentally than the 
neutron dripline, nuclear forces remain unexplored in 
medium-mass proton-rich nuclei. Because the proton 
dripline is closer to the line of stability, it has also been 
mapped out empirically by calculating the energies of 
proton-rich systems from known neutron-rich nuclei us- 
ing the isobaric multiplet mass equation (IMME) 0, [l(| 
or Coulomb displacement energies [111 ]. This suggests 
that proton-rich nuclei provide an important testing 
ground for nuclear forces including known Coulomb and 
isospin-symmetry-breaking effects. 

In this Letter, we present the first study of chiral 3N 
forces for proton-rich nuclei. The couplings of 3N forces 
are fit to few-nucleon systems only, and we provide pre- 
dictions for the ground-state energies (Figs.Q]and[3]) and 
spectra (Figs. [2] and 2]) along the chains of N = 8 and 
N = 20 isotones to the proton dripline. Our results show 
that chiral 3N forces provide key repulsive contributions 
in proton-rich nuclei, similar to the neutron-rich case. 
The precise location of the dripline in N = 8 depends 



on the calculation, with a potential two-proton-decay 
candidate 22 Si. For the N = 20 isotones, we robustly 
predict the dripline at 46 Fe and the two-proton emitter 
48 Ni 12|,|l3(. Furthermore, we find good agreement with 
experimental spectra of 18 Ne, 19 Na, 20 Mg and 42 Ti and 
provide predictions for the isotones where excited states 
have not been measured. 

We consider a shell model description of the N = 8 
and N — 20 isotones and determine the interactions 
among valence protons, on top of a 16 and 40 Ca core, 
based on nuclear forces from chiral effective field theory 
(EFT) [14]. At the NN level, we take the chiral N 3 LO 
potential of Ref. [H| and evolve to a low-momentum in- 
teraction Vjowfc with cutoff A = 2.0 fm -1 using renormal- 
ization group methods, which improve the many-body 
convergence [161 ]. Three-nucleon forces are included at 
the N 2 LO level. These consist of the long-range two-pion- 
exchange part, as well as one-pion-exchange and short- 
range contact terms [14j . The shorter-range 3N couplings 
Cd and ce are determined by fits to the 3 H binding en- 



ergy and the 4 He radius for A3N = A = 2.0 fm -1 [l7j |. 
without further adjustments in the many-body calcula- 
tions presented here. 

Excitations outside the valence space are included 
to third order in many-body perturbation theory 
(MBPT) [li| in a space of 13 major harmonic-oscillator 
shells. We have checked that the matrix elements are 
converged in terms of intermediate-state excitations. For 
the N — 8 isotones, we consider both the standard 
sd-shell and an extended sdf7/ 2 p 3 / 2 valence space with 
hcu = 13.53 MeV, and for N = 20, the pf and pfg 9 / 2 
spaces with Huj = 11.48 MeV. The extended valence 
spaces proved important in this framework for oxygen 
and calcium isotopes @-[l]. In addition to the NN- force 
contributions, we include the normal-ordered (with re- 
spect to the core) one- and two-body parts of 3N forces 
in 5 major shells 0, 0|- The normal-ordered parts domi- 
nate over the contributions from residual three-body in- 
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TABLE I. Empirical (emp) and calculated (MBPT in the 
standard/extended valence spaces) SPEs in MeV. 



teractions [5, The latter are expected to be weaker 
in normal Fermi systems due to phase-space limitations 
in the valence shell compared to the core :20] . 

For the valence proton single-particle energies (SPEs) 
in 17 F and 41 Sc, we solve the Dyson equation self- 
consistcntly including the contributions from NN and 3N 
forces in the same spaces and to the same order in MBPT. 
Our calculated SPEs are given in Table HI in comparison 
with empirical SPEs taken from the experimental spec- 
tra of 17 F and 41 Sc. The MBPT SPEs are similar to 
the empirical values, but the si/ 2 and py 2 orbitals are 
at higher energy in both spaces. This finding is similar 
to the calculated neutron SPEs in oxygen and calcium 
isotopes [1, 13] , which are more bound and differ due to 
Coulomb and isospin-symmetry-breaking interactions. 

All calculations based on NN forces are performed with 
the empirical SPEs in the standard sd- and pj-shells 
(NN forces only lead to poor SPEs), while those involv- 
ing 3N forces use MBPT SPEs in both standard and 
extended valence spaces. In this work, we focus on 3N 
forces, whose contributions are of the order of a few MeV, 
but an explicit inclusion of the continuum naturally be- 
comes important for weakly bound states and can lead 
to very interesting contributions, typically of several hun- 
dred keV [|, 21 1 . Therefore, we only show spectra to 22 Si, 
which is cither at the dripline or a two-proton emitter, 
and to the two-proton emitter Ni. Note that in the case 
of weakly bound or unbound states, additional attractive 
contributions from the continuum are expected 21 j . 

We first consider the ground-state energies of the N = 
8 isotones from 18 Ne to 24 S, which we compare with ex- 
periment when available. As data is limited, we also em- 
ploy the IMME 0] . This relates the energies in an isospin 
multiplet (of states with the same quantum numbers in 
different isobars A) by a quadratic dependence in isospin 
projection T z = (Z — N)/2, 



E(A, T, T z ) = a(A, T) + b(A, T) T z + c(A, T) T 2 . (1) 

The energies of proton-rich nuclei can thus be obtained 
from their — T z isobaric analogues via E(A,T,T Z ) — 
E(A, T, -T z )+2b(A/T)T z , using a standard fit of the em- 
pirical ^-coefficient @, b = (0.7068A 2 / 3 - 0.9133)MeV, 
with the atomic mass evaluation (AME2011) [22| for 
known neutron-rich nuclei. Moreover, we include for 
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FIG. 1. Ground-state energies of N = I 
to 16 0. Experimental energies (AME2011 
olations as open circles) and IMME values are shown. We 
compare NN-only results in the sd-shell to calculations based 
on NN+3N forces in both sd and sdf 7 / 2 p3/ 2 valence spaces 
with the consistently calculated SPEs of Table [I] 
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TABLE II. Experimental and calculated one- and two-proton 
separation energies S p and S2 P (in MeV) of N = 8 isotones. 
Where data is unavailable, IMME values are given in brackets. 



comparison the extrapolated values of AME2011, al- 
though the IMME is considered to be more accurate. 

Figure Q] shows the calculated ground-state energies, 
obtained from exact diagonalization in the valence spaces 
with NN-only and NN+3N forces, compared with the 
AME2011 experimental values and extrapolation, and 
with the IMME. As expected, the IMME reproduces well 
experimental data. It finds 22 Si to be bound, though only 
by lOkeV, with respect to 20 Mg. 

In the calculations based on NN forces only, we see 
a systematic overbinding throughout the isotone chain, 
which becomes more pronounced for larger mass num- 
ber. Three-nucleon forces provide key repulsive con- 
tributions to ground-state energies and good agreement 
with experiment is obtained in both valence spaces. The 
extended-space predictions become more bound beyond 
20 Mg, the last measured isotone. For the sd-shell, the 
proton dripline is predicted at 20 Mg, and in the extended- 
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FIG. 2. Excitation energies of N = 8 isotones calculated with 
NN+3N forces in the sdfy 2 p z / 2 valence space, compared with 
experimental data 23- 27] where available. 

space calculation, it is at 22 Si, which is only two-proton 
stable by 0.3 MeV. This makes an unambiguous predic- 
tion of the dripline difficult, and an experimental mea- 
surement of the 22 Si ground-state energy would present 
a decisive constraint for 3N forces. All calculations find 
a sharp decrease in binding energy past 22 Si, clearly in- 
dicating the proton dripline has been reached. 

A more detailed picture can be developed from the one- 
and two-proton separation energies given in Table HU S p 
and S^p are key quantities for determining two-proton 
emission candidates. The experimental (and IMME) sep- 
aration energies in most cases lie between the sd- and 
sdf 7 / 2 p 3 / 2 -space calculations. While the sd-shell energies 
are slightly closer to experiment for lighter isotones, the 
extended-space calculations agree best with the IMME 
beyond A — 20, approximately the same point, 21 0, at 
which the added valence-space orbitals become impor- 
tant in the oxygen isotopes in this framework [2]. 

Spectroscopic data in the N — 8 isotones exists to 
°Mg. In Fig. [5J we compare the experimental low-lying 
states in 18 Ne, 19 Na, and 20 Mg to those calculated with 
NN+3N forces in the sdf?/ 2 pa/ 2 valence space. Calcu- 
lations with 3N forces in the sd-shell give very similar 
spectra up to 20 Mg, while for 21 Al and 22 Si they are 
more compressed than in Fig. [2] In 18 Ne we find a very 
good agreement for the first excited 2 + and 4 + states. 
The ground state and first two excited states in 19 Na 
have been measured with tentative spin and parity as- 
signments 3, HI ■ The ordering of the first two states in 
our calculation disagrees with the tentative assignments, 
but the spacing between them is only 0.3 MeV. Further- 
more, while the calculated V 2+ energy is overpredicted, 
the Sy 2 orbital is unbound in our calculation, and contin- 
uum coupling is expected to reduce this energy. In 20 Mg, 
only information on the first excited state has been pub- 
lished [26 [jbut a second excited state has been measured 
recently 27[. While a tentative assignment of 4 + was 
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3. Ground-state energies of N 
40 r 



■- 20 isotones rela- 
tive to * u Ca. Experimental energies [28| (filled circles) and 
AME2011 extrapolations ^ (open circles), as well as IMME 
values are shown. We compare NN-only results in the p/-shell 
to calculations based on NN+3N forces in both pf and pfg9/ 2 
valence spaces with the SPEs of Tabic [I] 
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TABLE III. Experimental and calculated one- and two-proton 
separation energies S p and S2 P (in MeV) of N — 20 isotones. 
Where data is unavailable, IMME values are given in brackets. 
Direct measurements of S2 P in 48 Ni are from Refs. [13. f]~3j| . 



given to this state, we predict two close-lying states (2 + 
and 4 + ) at a very similar energy. In our predictions for 
21 A1 and 22 Si, we note the high 2 + state in 22 Si as a pos- 
sible indication of a subshell closure analogous to 22 O [2j. 

Next, we show in Fig. [3] the ground-state energies 
of the V = 20 isotones from 42 Ti to 48 Ni, where the 
IMME also reproduces well the limited experimental 
data . Calculations with NN forces already lead 
to a reasonable description of experiment, with ener- 
gies only modestly overbound (within 1 MeV) beyond 
45 Mn. When 3N forces are included, the additional re- 
pulsion systematically improves the agreement with data. 
The extended-space calculations agree very well with the 
IMME throughout the isotone chain, while the p/-shell 
results deviate for 47 Cr and 48 Ni. In all calculations the 
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FIG. 4. Excitation energies of TV = 20 isotones calculated 
with NN+3N forces in the p/g9/ 2 valence space, compared 
with experiment available for 42 Ti only [23l ]. 

proton dripline is robustly predicted at 46 Fe. 

The one- and two-proton separation energies are given 
in Table Mil The experimental and IMME values gen- 
erally fall within the NN+3N calculations in the pf 
and pfga/2 valence spaces. The difference in S p and 
S2 P between the two calculations only becomes larger 
than 0.5 MeV for 47 Cr and 48 Ni. This indicates that, in 
our framework, the <jy 2 orbital becomes relevant around 
A = 47 and provides extra binding, similar to the calcium 
isotopes 0, 3j- Indeed, our pfgs/ 2 result for S o v o f 48 Ni 
is only 0.3 MeV larger than recent experiment [13. Il3j|. 

Spectroscopic data is only available for 42 Ti in the 
N = 20 isotones. We show the predicted spectra based 
on NN+3N calculations in the pfgo/ 2 valence space in 
comparison with experiment in Fig. [4] The energies of 
the first 2 + , 4 + and 6 + are well reproduced. There are 
two observed states between 2f and 4^~ that do not ap- 
pear in our calculation. We attribute these to neutron 
(4p2h) excitations, expected around 40 Ca [2{|. For the 
remaining isotones, we show our predictions for the en- 
ergies of the first five excited states. Similar to 22 Si, we 
note the modestly high energy of the 2 + state in 48 Ni 
as a tentative closed subshell signature. The calculated 
spectra in the p/-shell are similar up to 44 Cr, and more 
compressed for the heavier isotones. 

In summary, we have presented the first study of chiral 
3N forces in proton-rich medium-mass nuclei. Our results 
for ground- and excited-state energies are in very good 
agreement with experiment, including the prediction of 
a recently discovered state in 20 Mg 27]. A future mea- 
surement of the ground-state energy of 22 Si will provide 
an especially sensitive test of chiral 3N forces. We make 
predictions for the unexplored spectra of the N = 8 and 
N = 20 isotones, which cannot be accessed by standard 
shell-model calculations with isospin symmetric single- 
particle energies. Our extended-space calculations for the 
ground-state energies are of similar quality as empirical 



IMME predictions, which is very promising for studying 
isospin symmetry breaking in medium-mass nuclei based 
on chiral EFT interactions. Our work presents a bridge 
to future studies, based on nuclear forces, of exotic nuclei 
with proton and neutron valence degrees of freedom. 
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